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Abstract 

The spherical reduction of the rational Calogero model (of type A n _ 1 and after remov¬ 
ing the center of mass) is considered as a maximally superintegrable quantum system, 
which describes a particle on the (n—2)-sphere subject to a very particular poten¬ 
tial. We present a detailed analysis of the simplest non-separable case, n=4, whose 
potential is singular at the edges of a spherical tetrahexahedron. A complete set of 
independent conserved charges and of Hamiltonian intertwiners is constructed, and 
their algebra is elucidated. They arise from the ring of polynomials in Dunkl-deformed 
angular momenta, by classifying the subspaces invariant and antiinvariant under all 
Weyl reflections, respectively. 


1 Introduction and summary 

The Calogero (or Calogero-Moser) model [1, 2, 3] 1 is the paradigmatical n-particle inte- 
grable system in one space dimension. Originally defined for the root system of A\ © A n _i, 
the Calogero model was quickly generalized for any finite Coxeter group of rank n [4], In 
particular the rational version is remarkable for its conformal properties and its maximal 
superintegrability [8, 9] (see also [10]). Moreover, when the coupling constant g is integral, 
the quantum model enjoys additional and algebraically independent conserved quantities, 
which make it ‘analytically integrable’ [11, 12]. In this case, intertwining operators relate 
the energy spectrum with that of the free theory (at g=0 or g— 1). This structure was 
worked out in the 90s by a Russian group [11, 12] using Darboux dressing and by Dutch 
mathematicians [13, 14] employing Dunkl operators [15, 16]. 

Maximal superintegrability for a quantum system of phase-space dimension 2 n has been 
characterized in the literature; for a recent account see [17] where it is termed ‘maximal 
degenerate integrability’. The Hamiltonian H is part of a rank-n abelian algebra of Liouvillc 
charges (commuting Hamiltonians), but this is embedded in a larger commutant of H, which 
is spanned by the Liouville charges plus n— 1 additional integrals of motion (Wojciechowski 
charges). The latter do not commute with one another but form a polynomial algebra [9] 
together with the Liouville charges. Denoting the commutant generators by C and distin¬ 
guishing between Liouville charges C' and Wojciechowski charges C" so that {C} = {C 1 , C"}, 
one has (see, e.g. [18]) 

[C',C'] = 0 , [C",C] = A(C, C") , [C",C"] = B(C’, C") , (1.1) 

where A and B are polynomial in C' and linear in C" , and we have suppressed the counting 
labels. This implies that each H eigenspace splits into simultaneous eigenspaces of the other 
n— 1 Liouville charges. Hence, we can find an H eigenbasis labelled by n quantum numbers 
(£ a ) with a = 1 , ,n. Furthermore, these simultaneous eigenspaces carry a representation 
of the algebra (1.1), which governs their degeneracies. 

Since the rational Calogero Hamiltonian is part of an SL(2,M) conformal algebra, it 
is natural to reduce the rational model to an integrable system on a sphere. Restricting 
ourselves to the customary Ai © A n _i model, one may employ the translational invariance 
and firstly separate the A\ part representing the center-of-mass motion. The relative-motion 
Hamiltonian may be interpreted as ruling a single particle in M n_1 , subject to a potential 
which is singular at the reflection hyperplanes and decays like r^ 2 with distance r from the 
origin. Secondly, this radial scaling suggests a further codimension-one reduction to S n ~ 2 , 
which defines the ‘angular (relative) rational Calogero model’. 2 This system describes a 
particle moving on S n ~ 2 in a collection of |n(n—1) Higgs-oscillator potentials [19, 20], each 
one centered at a positive root of A n _\ (and its antipode) and blowing up at the intersection 
of the reflection hyperplane with the unit sphere. These singular loci tessellate the (n— 1)- 
sphere into n\ hyperspherical simplices in which the particle is trapped. 

A number of classical features of angular rational Calogero model (center-of-mass-reduced 
or not) have been recently analyzed by an Armenian group [21, 22, 23, 24, 25, 26]. Intertwin¬ 
ing relations for the quantum model were investigated by M. Feigin already in [27], but more 

1 For reviews, see e.g. [4, 5, 6, 7]. 

2 It is investigated already in the appendices of [2] (see also [1, 3] for n= 3). 
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recently the energy spectra and eigenstates have been constructed [28], and the algebra of 
angular Dunkl operators has been studied [29]. These partial results leave open the following 
questions: 

• What is a complete set of algebraically independent conserved charges? 

• Can one construct a subset of Liouville charges (in mutual involution)? 

• What is a complete set of algebraically independent Hamiltonian intertwiners? 

• Can one find additional conserved quantities which get intertwined? 

• What kind of algebra is generated by the conserved charges and intertwiners? 

• Are these models maximally quantum super- or even analytically integrable? 

The answers are not easily obtained for the angular Calogero models because they are linked 
to nontrivial mathematical problems. For instance, the commutant of H is the ring of all 
Weyl-invariant polynomials in the angular Dunkl operators (for a definition see the following 
section), and it is not obvious how to construct a minimal set of algebraically independent 
generators. Moreover, to establish the superintegrability one should find the maximal abelian 
subalgebras (i.e. the Liouville charges) in this ring. 

Lacking a general strategy to attack these issues, we take an explicit look at the simplest 
examples, n =3 and n= 4. The first one, leading to a rank-one angular system, is quite trivial 
but well known, while the second one is still tractable and serves to illustrate most of the 
general concepts. Therefore, in this paper, we exemplarily analyze the rank-two quantum 
system which is obtained from the four-particle rational Calogero model, via translational 
reduction to M 3 and then radial reduction to S 2 . It describes a particle on the two-sphere 
trapped in one of 24 right isosceles spherical triangles which make up the spherical projection 
of a so-called tetrahexahedron. 

After presenting general (and mostly previous) results on the A n _i angular models in 
Section 2, the essentially trivial n=3 case, known as the Posclil-Teller model, is reworked in 
Section 3 as a warmup example. Section 4 introduces the n =4 model and recapitulates its 
energy spectrum and eigenstates [28], before our results concerning conserved charges and 
intertwiners are collected in Section 5. There, we establish that any conserved charge of the 
tetrahexahedric angular Calogero model is built from two independent operators, say J 4 and 
Jg (besides the Hamiltonian), and that this system features two independent Hamiltonian 
intertwiners, which we denote by M 3 and Mg. We provide explicit expressions for both and 
work out their (non-standard) intertwining relations with the elementary conserved charges. 
Since one cannot produce two commuting words with J 4 and Jg, the model enjoys only 
two charges in involution but, given three independent conserved quantities, it is maximally 
superintegrable. For integer value of the coupling, an additional independent conserved 
charge introduces a Z 2 grading (in the spirit of nonlinear supersymmetry without fermion-like 
degrees of freedom [30]) and renders the model analytically integrable. Various interesting 
algebraic relations are collected in the Appendix. 
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2 The angular (relative) rational Calogero model 


We parametrize the n-particle quantum phase space with coordinates x M and momenta p Ul 
subject to the canonical commutation relations (setting h — 1 ) 

[ x [L . p y ] = i 5 IJ U with n,u — 1,... ,n . ( 2 . 1 ) 


The center-of-mass coordinate and momentum, 


A' = and P = ^ p v 


[X,P] = i, 


( 2 . 2 ) 


decouple in the rational Calogero Hamiltonian 


H C ai 


+ E 


9(9 -!) 


= i-P 2 + H 


Li<v 


(x^—x 


v\2 2 n 


(2.3) 


and will be ignored in our starting Hamiltonian 


h = 


(2.4) 


which describes the relative motion of (A n _i) Calogero particles on the real line, or of a 
single quantum particle in R™ _1 subject to a potential singular at the A n _i Weyl chamber 
walls. The Hamiltonian is invariant under the reflection g —y 1 —g of the real coupling g, but 
higher conserved charges and intertwiners will not share this property. 

Since this Hamiltonian is scale covariant, it pays to introduce on M” -1 polar coordinates 
(r, 6) with a (relative) radius r and n—2 angles 6 together with their momenta p r and Lg, 

{ x v ) 2 = r 2 and \ ^(p^-Pu) 2 = p 2 r + ^L 2 + ■ (2-5) 

The (relative) angular momentum squared L 2 is best expressed by parametrizing the relative 
particle motion in M n_1 with coordinates y l where i — 1 , ..., n— 1 , 

n— 1 

7,2 = y ’ Pi = Py i , L ij = -i (y l Pj-y J Pi) , l 2 = ■ ( 2 - 6 ) 

i= 1 i<j 


Jacobi coordinates are a suitable choice for {y 1 }. Together with a dilatation generator D and 
a special conformal generator K, the Hamiltonian H forms an SL(2,R) commutator algebra, 
as is easily seen in the polar representation 


H 


12 , (n—2)(n—4) 

2 Pr "T" Sr 2 


+ ^Hr 


D = \{rp r + p r r) 


K = |r 2 . (2.7) 


Here, the angular (relative) Hamiltonian Hq is independent of the radial coordinate or mo¬ 
mentum, and it coincides (up to a constant shift) with the Casimir operator of the conformal 
SL(2,R), 


H n = \L 2 + U(6) = C — |(n—l)(n—5) with C = I< H + H I< - \D 2 . (2.8) 
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This defines our angular Calogero system (for A n _!), which describes a quantum particle on 
the sphere S n ~ 2 under the influence of the potential 


U0) 


r2 Y, 

H<u 


9 ( 9 - 1 ) 

(x^—x v ) 2 


r 2 

a£Tl + 


9 ( 9 - 1 ) 

(a ■ y) 2 


9(y9 2 ^ ^ cos 2 9, 


(2.9) 


where 1Z + is the set of positive roots for A n _i and 9 a denotes the angle (geodetic arc length on 
S n ~ 2 ) between the point 9 and the intersection of the ray in direction a with the unit sphere. 
Each of the \n(n— 1) contributions of cos -2 form is also known as a ‘Higgs oscillator’ [19, 20]. 
This potential is a very special one as it is tied to the root system. Its singular walls are of 
codimension one and tesselate the sphere into n\ hyperspherical simplices. 

For the spectral problem, let us pass to the position representation and encode momenta 
via differental operators, 


Pi 1 ^ -i di => p r -i(<9 r + %/) ( 2 . 10 ) 

so that 

H ^ -\(d 2 r + ^d r ) + o 

( 2 . 11 ) 

H n 1 y r 2 H + \(rd r + n— 3) rd r . 

It is customary to remove the first-order radial derivative term via a similarity transforma¬ 
tion, 

T(r, 9) = r 2 u(r,9) and H = r 2 H r 2 , (2.12) 

so that the Hamiltonian on R n_1 acts as follows, 

H'u ^ -\(d 2 r - (n ~ 2 $- 4) ) u + ±H n u . (2.13) 

As the potential is repulsive, the spectrum is positive semi-definite and continuous, 

H^ E)q = H'u E>q = E u E q with E e M > 0 , (2.14) 

where the wave function depends on further (suppressed) quantum numbers, but q 
parametrizes the angular Hamiltonian eigenvalues [28] (see also the appendices of [2]), 

H n v q = £ q v q with e q = lq (q + n — 3) and q G M > 0 • (2.15) 

For vanishing potential (g=0 or 1), = |L 2 , and q = £ E N 0 is the familiar total angu¬ 

lar momentum for a free particle on S n ~ 2 , with a degeneracy of (^™ 2 2 ) — (^n- 2 4 )- The 
interaction greatly reduces the degeneracy and shifts this quantum number [28], 

q = \n(n—l)g-\-£ and £ = 3£^ + 47 4 + ... + n£ n with £^ E No , (2.16) 

where the decomposition £ — > (£3, £4,..., £ n ) into additional quantum numbers yields the 
correct degeneracy, 


deg n (e g ) = p n (£) - p n (£-l) ~ p n (£-2) + Pn(£-S) , (2.17) 
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with p n (£) denoting the number of partitions of £ into integers not larger than n or, equiva¬ 
lently, the number of partitions of £ into at most n integers. Its generating function reads 

oo n 

Pn(t) := $>„(«)«' = (2.18) 

£=0 m= 1 

and the large-£ asymptotics is p n (£) = + 0(£ n ~ 2 ), implying deg n (e g ) ~ £ n ~ 2 for 

£ —y oo. Note that the degeneracy does not depend on g but only on £ (for any fixed n). Up 
to n— 5, fairly explicit formulae may be found [31]: 


deg 3 (f) 

deg 4 (£) 


deg 5 (£) 


0 for £ = 1,2 mod 3 

1 for £ = 0 mod 3 



0 for £=1,2, 5 mod 12 
1 for £ = else mod 12 


6£ 2 + 72£ - 89 
720 


{ 0 for £ = 2,22,26,46 mod 60 

2 for £ = 0,48 mod 60 

1 for £ = else mod 60 


(2.19) 


The wave functions may also be given explicitly. From 

(d 2 r -(q+%-l){q+%-2)r~ 2 + 2E)u E , q = 0 (2.20) 

we read off that 

U E ,q{r,0 ) = \fr Jq + ( n _ 3 )/ 2 (V 2 S r) v q (6) , (2.21) 

revealing a radial wave function of Bessel type. The angular wave function v q {6) depends 
on all additional quantum numbers (£ 3 ,£ 4 ,... ,£ n ) and is a bit harder to find [28]. It can be 
constructed in the following way, 



,(<?) 


(0) 


r n—3+q 


n 

/i=3 


n 

x> 

V=1 


v ) 


^9 r 3- n -n(n-l)g 


( 2 . 22 ) 


where 

A = n^-x") = H a y (2.23) 

g,<v a£ TZ+ 

is the Vandermonde of the original M n coordinates and T> u denotes the so-called Dunkl 
operator [15, 16], 


V u 


d u 


A 


1 


x v —x p 


&vp 


[V^V v ] = 0 , 


(2.24) 


with s up permuting the position and momenta of the uWi with those of the pth particle. 
The Dunkl operators play a central role in asserting the integrability of the Calogero model, 
because deforming d u —)• T> u essentially creates the interacting system from the free one. 
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Their Newton sums yield the Liouville charges of the Calogero system, and their mutual 
commutativity guarantees that of the charges. 

It is convenient to switch again to the relative coordinates {?/*}, in which the Dunkl- 
deformed momenta and angular momenta take the form (up to a factor of —i) 

= di - g ^ and C tJ = ~(y l Vj - y J Vi) , (2.25) 

a &n + a ' y 


respectively, where s a is the Weyl reflection on the hyperplane orthogonal to the root a. 
Polynomials in the Dunkl operators are not just differential operators but also act by Weyl 
reflections. Their restriction to functions totally symmetric under the Weyl group is called 
‘residue’ (denoted by ‘res’) and produces a purely differential operator. In this way, the 
Hamiltonian can be obtained from the homogeneous symmetric polynomial of degree two, 

« = = “IE®? =*• H = (2.26) 

i 

Analogously, the symmetric quadratic polynomial in the angular Dunkl operators Cij shifted 
by a certain pure-reflection term, 

W.n = ~\ ^ £ 2 , + \ S {S + n- 3) with S = g^2 Sa > (2.27) 

i<j a 

provides the angular Hamiltonian by taking the residue [29], 3 

H n = res (Hq) = |res(£ 2 ) + \ rt(n ~ 1) p ( n(n 2 ~ 1) p + n-3) = |res(£ 2 ) + e q (t= 0) . (2.28) 

We remark that the residue of the pure-reflection term simply produces the ground-state 
energy e q (£=0). Clearly, H and Hq may be considered as ‘Dunkl deformations’ of (—\ 
times) the Laplacian on M n_1 and S'”” 2 , respectively. 

To express the angular wave function v q {6) in terms of the relative coordinates, we must 
rewrite A in terms of {y 1 } and the Newton sums of the V v in (2.22) in terms of the Vi, 
which depends on our choice for {?/*}. Phrased differently, we need to identify the totally 
Weyl-symmetric homogeneous polynomials (7 /i ({'D,}) of orders fi = 3,4,..., n, to go inside 
the large brackets of ( 2 . 22 ), 

v q {9) = v^iO) ~ r n - 3+q ff[a2{Vi}Y^ A 9 . (2.29) 

V=3 ' 


From the form of (2.29) one can infer that v q (9) is a ratio of a degree-g homogeneous polyno¬ 


mial h q (y) in {y 1 } to the gth power of the radial coordinate r = [JA(y*) 2 ] • From (2.11), 

the r independence of v q {6) and (2.15) it follows that the polynomial h q (y) is annihilated 
by H, 


h[ g) := 


:= r* v 


(g) 


( 0 ) 


Hh[ 9) = 0 


(2.30) 


3 We note that C 2 = — J2icj ^ij i s a non-negative operator. 
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and hence may be viewed as a Dunkl-deformed harmonic polynomial on i? n ~ 1 . Another 
property visible from (2.29) is that h q (y) contains a factor of A 9 , which may be split off to 
define 

h {9) := A~ 9 h[ 9) =*► Hh {9) = 0 , (2.31) 

which then is in the kernel of the Hamiltonian ‘in the potential-free frame’, 

H = A~’H A» = + ^(9^)} . (2.32) 

The potential-free Dunkl-deformed harmonic polynomial }if } is homogeneous of degree ^ 
only, for any value of g. In particular, for the ground state one has 

h (9) = 1 => h {9) = A 9 , (2.33) 

and hence the full ground-state wave function is totally symmetric (antisymmetric) under 
particle permutations or Weyl reflections for even (odd) integer values of g. Since all other 
ingredients besides A 9 in ( 2 . 22 ) are completely symmetric, this symmetry property of the 
integer-g ground state extends to all excited states above it. Furthermore, the reflection 
symmetry g+ 1 O —g is broken since one tower of states is Weyl symmetric while the other 
one is antisymmetric. 

It is well known [13, 14] that there exists an intertwining operator M = M^ which 
establishes an isospectrality of and H^ 9+1 \ This differential operators of order |n(n—1) 
has the simple form 

M = res (M) with M = J J(V„-V v ) = JJ a V , (2.34) 

a£lZ + 


and using the Weyl antisymmetry of A4 it is straightforward to verify that (see e.g. [27]) 
[M,n\ = 0 => M {9) H {9) = H {9+]) M {9) and (2.35) 

so that 

M (s) : (g,£,q) (g+l , i, q+\n{n- 1)) (2.36) 

simultaneously acts on the radial wave function y/r J q+ ( n - 3)/2 and on the angular wave func¬ 
tion . Up to ro=4 the expressions for M ^ were worked out explicitly in [18]. 

However, the angular system Hq comes with its own intertwiners M s = Ms^ of some 
order seN, which then provides additional intertwiners for H. This is suggested by the form 
(2.16) of the quantum number q appearing in (2.15), which reveals a partial 4 isospectrality 
for H n , 

M { s 9) : (g,£,q) ^ (g+l, t-\n{n-l) , q) . (2.37) 

4 meaning that the intertwiners have a nontrivial kernel 
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Indeed, as shown in [27], one can more or less copy the strategy from R" 1 to S'" - 2 and 
obtain the lowest-order angular intertwiner, with s = \{n— l)(n—2), 


M n = Mi (n _ 1)(n _ 2) = res(AIo) with 

Mn = ^sgn(7r)7r( JJ {(x 1 -x 9 )(V 1 -V u ) ~ 

n£S n 2<[i<v<n (2.38) 

= x^ sgn( ^ 7r ) 7r ( n , 

n£S n 2 <fj,<v<n 


where the choice of x 1 as a reference is arbitrary and irrelevant. Because H is in the center 
of the algebra generated by {£p} and A4q is again antisymmetric under Weyl reflections, 
one argues, analogously to (2.35) and with the help of = H (a+1 \ that [27] 

[£«,«] = 0 =*■ and ~ *^> (n _ 1)/2 . 

(2.39) 

These properties hold not only for Mq, but actually any Weyl- antiinvariant function of 
{jCij} provides an angular intertwiner which obeys (2.39). By the same token, any 
Weyl -invariant function of {C t j} will yield, via its residue, a conserved angular quantity C t 
of some order t. We remark that the angular Dunkl operators do not commute with each 
other but form a subalgebra of a rational Cherednik algebra [29]. It is a deformation of the 
so{n— 1) Lie algebra generated by {Lij}. The general problem of identifying all Liouvillc 
charges for the angular system reduces to constructing a maximal abelian subalgebra in 
the algebra of Weyl-symmetric polynomials in {£y}, which is not an easy task. Likewise, 
identifying the minimal independent set of Weyl-antiinvariants and hence intertwiners is 
nontrivial. 

Now, since any such M s and Ct is scale invariant, i.e. r independent, it does not touch 
the radial wave function, and hence it is also true that 


as well as c[ s) 


and ~ 

(jig) 


To summarize, we have the connection 


C t (£ij) Weyl-invariant 
M. s (Cij) Weyl-antiinvariant 


C t = res (C t ) commutes with Hq , 

(2.41) 

M s = res(A4 s ) intertwines with . 


In general, the angular conserved charges C t will not be in involution with one another, but 
some combinations may be. The angular intertwiners M s are differential operators of order s 
and have, in contrast to M, a sizeable kernel, with 


dim ker n (£) = deg n (£) - deg n (£ - |n(n-l)) . 


(2.42) 


For small values of n, we can be explicit: 


dim ker 3 (£) 


dim ker^Q 


dim ker 5 (£) 



for £ = 0 
for £ > 0 


for £ = 0, 3 mod 4 
for £=1,2 mod 4 



for l — 0, 3,4, 5, 8, 9,11 mod 12 
for £= 1,2,6,7,10 mod 12 


(2.43) 


With the help of the adjoint intertwiners, which reverse the direction of M resp. M s , 

M (») t = M (-») => and ~ o“,(„-i )/2 , 

(2.44) 

it is trivial to construct special conserved charges, 

4? := and R { 2 ~ a) = M^M^~ 9) , (2.45) 

which provide particular cases of some Ct for t=2s. In the full Calogero model, R = M 
turned out to be a specific polynomial in the Liouvillc charges, which was evaluated by 
infinitely separating the particles thus making contact with the free case. This move is not 
possible for R 2s due to the compactness of S n ~ 2 . Neither is it true that R^ 9 ' 1 = i? 2 s +1 \ and 
so these charges intertwine as follows, 

= R£ 3) M<‘1 = (fl'f 11 - (l+2 9 )p<f ■>) Mj»> . (2.46) 

and similarly in the opposite direction. Here, p 2s is a particular expression in other conserved 
charges. 

Quite generally, we cannot expect conserved charges other than the angular Hamiltonian 
to be invariant under g —> 1 —g, and in fact one finds that 

C[~ 9) ± C[ 9+1) and M { s 9) C { t 9) ± c[ 9+1) . (2.47) 


Therefore, higher conserved quantities Ct intertwine in a more complicated fashion, 

, (2.48) 

S’,t' 

where T s s 'f (g) are certain polynomials in g of order s+t—s'—t'—l at the most. The sum runs 
over all linearly independent conserved quantities including 

C {9) = -res(£ 2 ) and C { 0 9) := 1 . (2.49) 

When the coupling g is an integer, repreated intertwining relates all quantities with their 
analogs in the free theory (g =0 or 1), which allows one to generate analytic expressions 
for all wave functions. Moreover, in this case additional conserved charges make the model 
‘analytically integrable’ and produce a ‘supercomplete’ graded ring of commuting differential 
operators [11, 12, 18]. 
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3 Warmup: the hexagonal or Poschl-Teller model 


Let us work out and test the general assertions of the previous section on the simplest 
nontrivial case, which occurs for n = 3. The three-particle (or A 2 ) rational Calogero model 
is long known to be completely separable, the ensuing angular (relative) model being the 
famed Poschl-Teller model of a particle on a circle with a cos -2 (3(f>) potential. 

We start by introducing Jacobi coordinates (y 1 , y 2 ) for the relative motion (and the center 
of mass X ), 


x 1 = X + ^y l + ^y 2 , d x i = ld x + ±d y i + ±d y 2 


= x -^y 1 + 7Ey 

x3 = X ~Ts y2 ’ d * 3 = y )x ~7t> d 

and then polar (r, </>) and complex (w,w) coordinates on M 2 , 
y 1 = r cos (j) and y 2 = r sin 
After the standard similarity transformation, the (reduced) Hamiltonian 

+ 


1 „,2 


d x 2 = \d x - -7s d v i + -4 d v 


V2 y 1 Ve y 


(3.1) 


w := y 1 + \y 2 = re 1 ^ 


(3.2) 


H ~ -Wi + dl) + 9(9-1) + J^TZ 

= -m+di) + 


(y 1 ~V3y 2 ) 2 (y'+VZy 2 ) 2 

9^(^-l)((2/ 1 ) 2 +(?/ 2 ) 2 ) 


(3.3) 


2 {y 1 ) 2 (y 1 —V3y 2 ) 2 {y 1 +V3y 2 ) 2 


and its angular cousin take the following form as differential operators on M 2 and S 1 , respec¬ 
tively, 

H'u ^ -i(9 2 + i ^)u + H n u with H n = —dfy + U, (3.4) 
where the angular potential receives contributions from the three positive roots of A 2 , 


U(4>) = 2 (cos 2 (</>) +cos 2 ((/>+^) + cos 2 (0—^)) = |#(g-l) cos 2 (30) . (3.5) 

The poles separate the circle into six disjoint arcs, pertaining to the hexagonal symmetry of 
the A 2 root system. We display the angular Hamiltonian in complex coordinates as well, 


Ho = \(wd w -wdyj) 2 + g(g- 1 ) 


18 (ww) 3 
[w 3 + w 3 ) 2 


(3.6) 


Specializing (2.15)-(2.17) to n= 3, one finds that the angular momentum t = 3£ 3 is a 
multiple of three and 


£ - In 2 

tq — yi 


with q = 3 g + i = 3(g + f? 3 ) and deg(e 9 ) = 1 . (3.7) 

The y/r factor in (2.21) gets cancelled in the similarity transformation (2.12), giving 

'b E , q (r,4 >) = J q (V 2 Er) v q {4>) , (3.8) 
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and the angular wave function reads 


v q (<j>) = v { e 9 \(i)) ~ r q A 9 r ~ 69 ~ r q (V 3 W - V% ) h A 9 r~ e ° (3.9) 

u=l 

with 

A ~ (x 1 —x 2 )(x 2 —x 3 )(x 3 —x 1 ) ~ r 3 cos(30) ~ w 3 + w 3 , (3.10) 

where we have rewritten in complex coordinates. 5 Owing to the S 3 Weyl group of A 2 , we 
take profit from the complex cubic roots of unity to cast the holomorphic Dunkl operator 
D w = -;(Vi — YD 2 ) into the simple form 


D w = d w - g{ —s 0 H - 3 S+ + -— 7 — : s_j with p = e 2m/3 , 

l w + w pw + pw pw + pw ) 

where the three basic Weyl reflections act as follows in the complex plane, 

So : w 1 —> —w , s_|_ : w 1 —^ —pw , s_ : w (->• —pw , 

or on the polar angle, 

s 0 : 0 h- > 7r — (j) , s + : 0 h- > — | — (j) , s_ : </> 1 —> | — <t> ■ 
According to (2.34), the full intertwiner is a third-order operator , 6 
M ~ (D 1 -V 2 )(V 2 -V 3 )(V3-V 1 ) ~ Dl + Vl, 

whose evaluation eventually produces [18] 


M 


d\->d‘yidu — -ppr d‘>:\d> I — <9l2 — Z^rdiodo 


I 4g 2 O I 4g 2 p, , 4g 2 s(ff-l) ^ fl(s-l) O 

+ ZTZZ23 O 31 + C'12 + Z3TZT1 <^23-7VTT\2 <A2 — f„S3l2 ^23- 7Z5TV2 @31 


(S 12 ) 


(S 23 ) 


(i 31 ) 


- ’/C'.','. + 2j(a-l)(a+2)(pip + pip + pip) . 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


with the abbreviated notation x MI/ = x^—x v and d^ u = d x »—d x v. In complex form, it reads 


M 


d 3 + d 3 - 

W ' W 


6 g 


w 3 +w 3 


(w 2 dl + ww d w du, + w 2 dl) + 


6g(3ff-l) 

w 3 +w 3 


(wd™ + wdu,) 


1) f 4 rj _ 4 rj ^ 3g(9g+14) 27g(g-\-2) ^ 6 _ 6 0 3 -.-p 

■ (ro + ro + - (ro +«; + 2 S ru> m ) , 


(3.16) 


(tc 3 +w 3 ) 2 


rc 3 +w 3 (tc 3 +tc 3 ) 3 


and in the polar representation the intertwiner becomes 


, 3(9g+l) 

I r 2 


M ~ cos 3 4> j<9 3 — 3(3g +^ d 2 

— sin 3 (j) Mfi | 3 d 2 — ^d r + ^ + | , 


dr-^ + ^H n d r - b -^H n 


(3.17) 


5 We refrain from displaying v q ((f>) in Jacobi or polar coordinates as the expressions are not enlightening. 

6 We remark that 2?^ + TA changes sign under any Weyl reflection while 2?^ — VA is totally symmetric. 
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where the angular intertwiner Mq (to be given shortly) has been used. 

Reducing to the angular subsystem, i.e. the Poschl-Tcller model, the formula (2.38) yields 


M.q = M-i ~ x 1 V 2 — x 2 V 1 + x 2 V 3 — x 3 V 2 + x 3 T>i — x 1 V 3 ~ i {wV w — wV^,) , (3.18) 

which simplifies to 

, „ o cx 1 +x 2 —2x 3 x 2j rx 3 —2x l x 3 +x 1 —2x 2 i , 

Mi ~ x 1 d 2 3+x 2 d 31 +x 3 d 12 + g { - — 2 -s 12 d- 35 -s 2 3 H--y- s 31 \ (3.19) 


x 


x 


23 


X 


31 


or to 


M± ~ i (wd w -wdu,) - i g 
Taking the residue, we arrive at 


w — w pw — pw pw — pw 

-v So H-+ —-— S_ 

w + w pw + pw pw + pw 


(3.20) 


, „ ~ .. „ (x lj rx 2 —2x 3 )(x 2j rx 3 —2x 1 ){x 3 -\-x l —2x 2 ) 

M n = M ! ~ x 1 ^ + x 2 0 3 i + x 3 d 12 - g - A 19 9 , ,, - -- 


~12 -23 ™31 

w IV IV 


(3.21) 


or, in complex and polar parametrizations, 


w 3 — w 3 


Mi ~ i (wdyj — wdyj') — 3i g — — — = + 3g tan30 = cos 9 (30) <9^ cos 9 (30) . (3.22) 


The penultimate expression is the familiar first-order operator which enjoys the following 
remarkable properties: 

• it adds zeros at 0 = 0, ^ and poles at 0 = n, ^ 

• it intertwines between symmetric and antisymmetric Weyl group singlets 

• acting ^ 3+1 times on v'^ gives zero 

• acting with m\ 1 ^ = M\ on n® yields non-normalizable wave functions v ^ _1 ' ) 

• it is a Vandermonde dressing of the angular momentum 

• it is the unique angular intertwiner 

For g— 0, the expression (3.9) is readily evaluated (with r° —» lnr and £ = 3£ 3 ), 


,(o) 


r l — d%Y z In ww ~ r e (w e ±w e ) 


cos £0 for £ even 
sm £6 for £ odd 


(3.23) 


which are all Weyl symmetric. Acting repeatedly with Mi on these wave functions, one al¬ 
ternates between Weyl symmetric and Weyl antisymmetric and generates the wave functions 
at the same value of £ for any integer coupling. The low-lying states for g = 0,1, 2, 3 are 
tabulated in Table 1. The only exception is the one-dimensional kernel of Mi, consisting of 
the ground state 

v Y' > ~ r~ 3g A 9 ~ cos 9 30. (3.24) 

Finally, the additional conserved quantity R 2 manufactured with Mq and its adjoint, yields 
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q — 3g + 3£ 3 
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£ = 0 g = 1 <7 = 2 

Figure 1: Poschl-Teller model: spectrum and intertwiner actions, 
nothing new since 

R 2 ] = (dc/, — 3g tan 3<j>) (d# + 3g tan 30 ) = — 2H^+9g 2 = —res(£ 2 ) = —C 2 \ (3.25) 
which provides p 2 l+l ' 1 = 9. We read off the intertwining relation 

m[ 9) R { 2 9) = ( R ( 2 9+l) - 9(1+2 g)) m[ 9) . (3.26) 

For generic values of g , there can be no further independent conserved quantity. How¬ 
ever, when g is integer, one additional algebraically independent conserved charge can be 
constructed, 

Q ia) = m[ 9 ~ 1) m[ 9 ~ 2) ■ ■ ■ M^~ 9) ==► Q {9) H {9) = = H^Q (a) , (3.27) 

which is a Weyl-antiinvariant differential operator of order 2g—1 that squares to 

(0 (s) ) 2 = II (-2< + 9f) . (3.28) 

1 = 1-9 


£ = V 
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4 Tetrahexahedric model: the spectrum 


The main purpose of this paper is to work out the angular (relative) rational 4-particle (or 
A 3 ) Calogero model as the simplest non-separable case, in order to demonstrate the viability 
of the general observations of Section 2. We take advantage of the isometry A 3 ~ D 3 and 
employ -D 3 -type relative coordinates . 7 

To construct the model and find its spectrum, we introduce not Jacobi but Walsh- 
Hadamard coordinates {y 1 } and their derivatives together with the center of mass X via 

x 1 = A + |(+t/ 1 + y 2 + y 3 ) , d x i = \d x + \{+d y i + d y 2 + d y 3 ) 

x 2 = X + \(+y l - y 2 - y 3 ) , d x 2 = \d x + \{+d y i - d y 2 - d y3 ) 

x 3 = X + \{-y x + y 2 - y 3 ) , d x3 = \d x + \{~d y i + d y 2 - d y3 ) 

x 4 = X + \(-y 4 - y 2 + y 3 ) , d x 4 = \d x + \{-d y3 - d y 2 + d y3 ) 

The polar (or spherical) parametrization of R 3 defines the radius r and two angles 9 and 4 >, 
y 1 = r sin# cos</> , y 2 = r sin 9 sin</> , y 3 = r cos 9 . (4.2) 

To avoid index cluttering, we redenote the Euclidean R 3 coordinates in a traditional fashion, 


(4.1) 

1 


(y 1 , y 2 , y 3 ) =■ (x,y,z). 

The angular momentum components read (up to a factor of i) 

L x ( yd z y ) , L y (zd x x<9 2 ) , L z (x<9^ y&x') > 

so that (minus) the Laplacian on S 2 takes the form 

= -(Ll + Ll + Ll) = -^sin 9d e --^d 2 , 

and the (reduced) Hamiltonian becomes 


H — ~T 2 {dl + d 2 + d'l) + 2 < 7 ( 17 —1)^- 


x 2 + y 2 y 2 + z 2 

~ t - / O O \ O 4” 


z 2 + x 2 


. (x 2 — y 2 ) 2 (y 2 — z 2 ) 2 (z 2 — x 2 ) 2 

After the standard similarity transformation, we obtain 


H' u 1 —>• — ^d 2 u + -^Hqu and Hq = ^L 2 + U 


(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 


with 

TT/ . . . f 1 cos 2 6 + sin 2 9 cos 2 <f> cos 2 9 + sin 2 9 sin 2 <f) 1 

U(0 , 0) — 2g(g—l)< — — —— + - — : y- 2 ±\2 4" 7 T7 ; TT • TXA 1 ' 

(sm 9 cos z 20 (cos z 9 — sm 9cos z (p) z (cos z 9 — sm 9 sin <p) ) 

(4.8) 

This has been named ‘cuboctahedric Higgs oscillator’ potential [21] since it is the superpo¬ 
sition of six Higgs oscillators placed at the positive roots of A 3 ~ D 3 , and those (together 

7 See also [32] for an application of the same idea. 
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with their antipodes) form the vertices of a cuboctahedron [33]. However, a better name 
is ‘tetrahexahedrah potential, because the intersections of the Weyl chamber walls with the 
two-sphere form six great circles, which trace the 36 edges of a tetrahexahedron (or tetrakis 
hexahedron) [34]. Its 24 identical faces are right isosceles triangles. The spherical projec¬ 
tion tessellates the two-sphere into 24 spherically triangular domains in which the particle 
is trapped. This sequence of geometrical configurations, starting from the A 3 root system 
(cuboctahedron) and ending at the walls of the angular potential (tetrahexahedron) is dis¬ 
played in Figure 2, and color and contour plots of the potential are shown in Figure 3. 



Figure 2: From the cuboctahedron to the tetrahexahedron. 




Figure 3: The tetrahexahedric potential (its log log log) on the two-sphere. 
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The angular energy eigenvalues become 

£ q ~ \9 (?+l) with q = 6 g + £ = 6 g + 34+44 . (4.9) 

The interaction reduces the 24b 1 dimensional free eigenspace (spanned by the spherical 
harmonics Y £rn ) to the subspace invariant under the S 4 Weyl group action generated e.g. by 

s x+y : (x, y, z) {-y, -x, +z) , s x - y : (x, y, z) (+y, +x, +z) , 

Sy+z : (x,y,z)t-> (+x,-z,-y) , s y - z : (x,y,z) (+x,+z,+y) , (4.10) 

s z +x ■ (x, y, z) (-2, +3, ~x) , s z - x : (x, y, z) 1-+ (+2, +y, +x) . 

This reduces the degeneracy to deg(e: 9 ) = deg 4 (£) of (2.19). The corresponding wave func¬ 
tions 

^E, q (r, 0, 4>) = j q (Y 2 Er)v q (e,(f)) (4.11) 

have a spherical Bessel function as their radial part and involve two specific Weyl-symmetric 
homogeneous Dunkl polynomials, of degree three and four, 

v q (9,<l>) = ^ 9) (M) ~ r q+1 3 A 9 r~ l ~ 129 


fl=l l/=l 

+3 


(4.12) 


~ r q+1 ( V x V y V z y (V A x +V A y +Vt) 4 A 9 r 1 - 129 , 
as well as the coordinate Vandermonde, 

A = \{m-W = (^ 2 - y 2 )(x 2 - z 2 )(y 2 - z 2 ) 

i<j 

~ r 6 sin 2 9 cos 4 9 cos 2 20 (tan 2 9 cos 2 0 — 1) (tan 2 9 sin 2 0 — 1) 
Here, the three Dunkl operators on M 3 are given by 


(4.13) 


V x = d x 


^y — dy 


v x = a 


9 „ 9 

■ s x+y 

x+y x—y 


&x—y 


9 


9 


9 


^x+y 

y+x y—x 


s x—y 


x+z 

9 


$x+z 


9 


x—z 

9 


®z+x 


<5?— 7 


1 °y+z 

y+z y-z 

9 9 


1 s y+z 

z+y z-y 


Sy—z 5 


s y-z 


(4.14) 


z+x z—x 

The angular wave functions u 9 (0,0) in (4.12) are known as ‘tetrahedral harmonics’ in theo¬ 
retical chemistry [35, 36]. 8 The low-lying states for ^7 = 0,1,2 are tabulated in Table 2. 

For completeness, we here display some quantities in the potential-free frame. The 
potential-free Hamiltonian 

H = -1(^ + 3’ + ® - + + + (4.15) 

annihilates the corresponding Dunkl-deformed harmonics 

h ( ^ (4 ~ r .i 2 s+i+ 6 t3+ 8 4 A —a (p x V y V z ) (:i (V 4 r +Vl+VT) U A a r'-' 2g , (4.16) 

so that the complete wave function may be written as 

^E/ 3 /S x ’y^ z ) = j g (V^r)r _ 5 A ff ^ 4 (x, 2 /,z) with 3 = 63 + 34 + 44 . (4.17) 


s The chemistry literature seems to be unaware of the straightforward construction scheme (4.12). 
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5 Tetrahexahedric model: intertwiners & integrability 


According to (2.34), the four-particle intertwiner is 

M ~ res((l>2 - Vl)(V 2 y - V 2 )(V 2 - V\)) , (5.1) 

and it acts on the wave functions as 

M {9) j q (V 2 Er)v < f\d,(j)) ~ j q+6 (V 2 E r) v^ +V> {9, <j>) (5-2) 

The rather complicated explicit form of this six-order partial differential operator has been 
constructed in [18]. Let us turn our attention to the angular intertwiners M s . Their building 
blocks are the Dunkl-deformed angular momenta 


C z = C yz =-(yV z -zV y ) , Cy = C z% = -(zV x -xV z ) , C z = C xy = -(xV v -yV x ) 

(5.3) 

or their complex combinations C± = C x ± i C y . Employing (4.4) and (4.14), they read 

£>X — L x + g{-^S X -y — -^S x+ y — : ^Sy- Z + y^Sy +Z + J^;S Z - X + -^-S z+X } , 

Cy = Ly + g{^^Sy- Z ~ ^^Sy +Z ~ J^S Z - X + J^S Z+X + -^S X -y + ^^S X+ y} , (5.4) 


£ z - L z + g{jz^ s z-x — J+^ S z+x — + ZTTT.Sx+'u + TT^Sy-z + 


1 


z-\-x^ z ~^~ x x—y^ x V 1 x-\-y X ~^~V y—z V z z-\-y V~^~ z J 

and transform under the Weyl reflections in the following fashion, 

$x-\-y • ^^xi ^yt ^z) 1 ^ (~^~^yi ~^~^x•> ^z) •> $x—y • (^#5 ^yi ^z) 1 ^ ^yi ^xi ^z) 1 

Sy-\-z • ^yi ^z) 1 ^ ( &x 1 ~\~^zi 1 $y—z • i^^x 5 ^yi ^z) 1 ^ ( ^xi ^zt ^y) ? (5.5) 

$z-\-x • {C X , Cy, C Z ) H > {-\-C z , — Cy 1 H - ^x) ? $z—x • x^xi ^yi ^z) 1 ^ ( ^zi ^yi ^x) • 

For convenience, we also note their transformation under certain even Weyl group elements, 

Sx+y^x—y • (^tc? ^yi ^z) 1 ^ ( ^xi ^yi^^z) 5 

Sy-\-zSy— z • (^£5 £z) ^ { + C X , — Cy, — C z ) , ( 5 . 6 ) 

Sz-\-xSz—x • (C X ,Cy,C Z ) I y ( — C x , +C y , — £ z ) . 

The construction (2.38) then yields a simple cubic polynomial in the angular Dunkl 
operators, 

M.Q. = AI3 = ^C X CyC Z -\- C X C Z Cy + CyC Z C X + CyC X C Z -\- C Z C X C y + C Z C y C X ^ , ( 5 . 7 ) 

which indeed is an anti-invariant under the Weyl group. It remains to evaluate the residue 
Mq = res(A4n) for arriving at 

M 3 — \(L x L y L z + L x L z L y ) — 2 g x $ 9 y 2 ( L x L y + L y L x ) 

+ 1G 9 2 - 4g(g-l) {y 2 t z if +2 g{-^s ~ ?rbO V zL x 

+ 2g , (fi'-l)(^+2) x 2 ^^2)2 +^(( y i z )3 - (^+2)3) + fi 1 ( 2 fi ,2 +8fi f_ l) §2^2 


+ 2 g 2 (8+90) (^23^)3^2^21 - |ff 3 ( B 2 - ya )(?_ 6 $)( y 2 - z 2 ) + cyclic permutations . 

(5.8) 
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Inserting the free angular momentum differential operators (4.4), we finally obtain 


44.3 U zd zxx yz d XX y + 2 (y z~')d xx + 4 g 2 _ z2 (yzd xx zxd X y 4 - x dy Z xyd zx j 


2 V 2 Z 2 ( __|_14a_ 

y Z \{x 2 -y 2 ){z 2 -x 2 ) + ( z 2 -x 2 )(y 2 -z 2 


-z 2 ) ( x 2 -y 2 ) 


2 9 -l , 2g—1 ^ 

-Z "I" ,_ 2 _ 2X J J 


(; z 2 -x 2 Y 


1x 2 y 2 


+ 


2x 2 z 2 


{z 2 —x 2 ) 2 (x 2 —y 2 ) 2 

„2 


2+ 


2z 2 


- 2 




xcb 


+ 2tf(0-l)(0+2) x 2 J 2 J^ 2) 2 + 2 + 9 ( 20 2 + 80 -l) jjszf 


(y 2 -z 2 ) 

„2^o 1 n„3 x 2 y 2 z 2 2 „3 a: 6 -H/ 6 +z 6 


+ 2 9 (8+99)^^^^^^ - 59 + c y dic p ermutations - 

(5.9) 

This is not the only elementary intertwiner however. At order six, one finds the next 
Weyl-antiinvariant function of {£ x , C y . £,}, namely 


Its residue then produces a second intertwiner, 

M, = {L‘,Lp+4g(-+ ; - J S ? ){^,Lp+4 9 (^ + lfc|l){Lj,Lp 


\y+z 
-2 +23 


+ 0(L 4 , L s x L y , L 2 x Ly) ± permutations . 


(5.10) 


(5.11) 


We have checked that all higher-order intertwiners, generated from further Weyl- 
antiinvariants, can be decomposed into linear combinations of M 3 and Mg with Weyl- 
synmretric functions of {C x , C y , C z } as coefficients. In this sense, M 3 and M 6 algebraically 
generate the set of all intertwiners for Hq,. As a curiosity, we report a remarkable (non-local) 
example from this set, 


M 


non-local 


= sin £ + sinh £ + + sin £_ sinh £_ = ^(— l) k 


£ 


4k+2 


+ £' 


4fc+2 


fc =0 

— i£ 10 +£ 10 l 

90 V + — J 1 113400V 

In all cases, the action on the angular states is as expected, 


= (4 + £ 2 _)-^(£ 6 + +£ 6 _) + 


(2fc+l)!(4fc+l)!! 
_ l _('£ 14 +£ 14 i 

681080400 V^-t- ) 




+ . . . . 
(5.12) 


(5.13) 


and the size of the kernels dim ker M s = dim ker 4 (4) is independent of g as was given in (2.43). 

What about additional conserved charges in our tetrahexahedral angular Calogero sys¬ 
tem? The S4 Weyl group action (5.5) on the Dunkl angular momenta is sufficiently simple 
to characterize the algebra of Weyl-symmetric polynomials. It is generated by 

J k := res(£*+ £* + £*) for k = 0,2,4, 6, (5.14) 

and its center is formed by 

Jo — C 0 — 1 and J 2 = -C 2 = -2 H n + 60(60+1) , (5.15) 
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Figure 4: Tetrahexahedric model: spectrum and intertwiner actions. 


so [+ 4 , Hq] = [+ 6 , Hq] = 0 but [+ 6 , + 4 ] 7 ^ 0. Any word formed from + 4 and + 6 will commute 
with Hq. The commutator [+ 6 , + 4 ] cannot be expressed in terms of lower charges but is 
related to a Weyl-invariant product of the two intertwiners, 

[J 6 ,J 4 ] = 12M|M 6 + 24(3+4#)+ 6 + 2 — 12(3+2 #)+ 4 + 4 — 48(l+2#)+ 4 + 2 + 2 
+ 12(1+2#) J 2 J 2 J 2 J 2 + 24(97+113#-96# 2 -144# 3 )+ 6 
+ 12(266+349#-304# 2 -528# 3 ) +4+2 + 12(76+101#-80# 2 -144# 3 ) + 2 + 2 + 2 
— 48(134+275#—1015# 2 —1596# 3 +1008# 4 +1728 # 5 )+4 
+ 24(146+379#—690# 2 —1240# 3 +672# 4 +1152# 5 )+ 2 + 2 
+ 96(10—39#—775# 2 —2232# 3 —1584# 4 )+ 2 • 

Higher-degree invariants are not independent, e.g. 


6+ 8 — 8+ 6 + 2 + 3 + 4+4 — 6 + 4 + 2+2 + + 2 + 2 + 2+2 

— 12(8+5#+12 # 2 )+6 + 4(34+23#+30# 2 )+4+ 2 — 8(5+3#+3# 2 )+ 2 +2+2 
+ 24(13+15#—102# 2 —72 # 3 )+4 - 4(43+70#-252# 2 -144# 3 )+ 2 + 2 

— 48(1+3#)(1+4#)(1—12#)+ 2 . 
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It is instructive to work out the complete set of elementary intertwinings: 


(9) 


{J>r } - 6(7+12 g)) M 3 
(j[ 9+1) - 4(11+12 g)4 9+1) + 48(26+73c/+48^ 2 )) M^ a) + 2 M f 


M^ 9 ) 4 9) = 

M^ 9) J {9) = 

M^ 9) 4 9) = (4 9+1) - (35+36#) j[ 9+1) -3(7+4g) J^ +1) J^ +1) 
+ 2(1111+2668^+1392c/ 2 )J^ +1) 

+ 96(A57+1933g+2717g 2 +1368g 3 +UAg 4 ))M^ 9) 

+ ( 34 9+1) - (115+200c/+48c/ 2 ))M 6 (9) , 

(■ 4 9+1) - 6(7+12^)) Mg 9) , 


(9) 


m 6 (9) + 2 (9) 


(5.18) 

(5.19) 

(5.20) 

(5.21) 


m 6 (9) +| 9) 


(—144Jg 9+1) + IUJ {9+1) 4 9+1) 


32 + 


(9+1) 7(9+1) 7(9+1) 

U n U o 


(5.22) 


7(9) 


- A8{53+75g+36g 2 )J { 4 9+1) + 16(77+99c/+36^ 2 ) +j 9+1) J { 2 9+1) 

- 96(261+1151g+1320g 2 +A32g 3 )J^ 9+1) 

+ 2304(159+1352(7+3672g 2 +4269g 3 +2232c/ 4 +432g 5 )) M 3 (9) 

+ (+j 9+1) - 4(29+36 g)4 9+1) + 72(27+139^+160/+48c? 3 )) M 6 ( 

with jjf 1 displayed in the Appendix. They have the general structure predicted in (2.48). 

Any maximally abelian subalgebra is generated by the center plus one arbitrary further 
element, hence only one Liouvillc charge besides the angular Hamiltonian can exist. There 
is no preferred choice among the words in + 4 and + 6 , but distinguished combinations are 


-Rfi = 


R 12 = 


m\m.3 


MjMg 


— \j& ~ \JiJ2 + I+2+2+2 — |( 11 + 165 i — ^g 2 )Ja 

+ +(13+245— 485 2 )J 2 +2 + 4(1+35)(1+45)(1—125)J 2 , 


12 + 6+6 + 12 {+ 6 , + r }+2 - y + 6 ' 72 ' 72+2 + 2 + 4 + 4+4 — 14 + 4 + 4+ 2 + 2 

lower-order terms , 


(5.23) 


6+4+2+2+2+2 — ~+2+2+2+2+2+2 + 


where the full expression can be found in the Appendix. The corresponding intertwining 
relations involve a single intertwiner only and read 

M^ 9) R {9) = [R {9+1) - (1+2 g)p { 6 9+1) ) M { 3 9) and 

M^R^ = (R[ 9 2 + l) - (l+2g)p[ 9 2 + 1} ) M< 9) with 

4 9+1) = -4+i 9+1) + 4+^ 9+1) + 2 (9+1) - 2(69+1285+485 2 )+^ +1) + 3{3+Ag)(7+l2g){ll+l2g) 

12+(9+i)+( 9+1)+(9+1) _|_ 294jl9+i) j-(9+i) 


36 +6 


1 ( 9+!) 
4 r 12 


= —54{+g 9+1) , +j 5 ' +1) } 


(9+1) 7(9+1) 7(9+1) 7(9+1) 

U Q L/O 


172 Jr ’ J, 


62+ 2 s+1) + 2 9+1) + 2 9+1) + 2 9+1) + 2 9+1) + 0(5) • 

(5.24) 

Similarly, the combinations M 3 M 6 and MgM 3 intertwine with a single intertwiner on either 
side of the relation (albeit not the same one). 
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Because the tetrahexahedral system is just two-dimensional, the existence of two con¬ 
served charges in addition to the Hamiltonian establishes its superintegrability. Beyond this, 
one can construct, for integral values of g, further conserved quantities, 

q( 9 ) = =*► Q {9) H[= Q {9 ) h£~ 9) = , (5.25) 

where M* stands for either M 3 or Mq, like it was achieved for the full Calogero model in [18]. 
These introduce a 7L 2 grading, with Q being odd and all previous charges being even. Clearly, 
the Q charges are algebraically independent of the even ones, but we expect that any two 
of them are related by some combination of even charges, and so they represent only a 
single additional algebraically independent conserved quantity. The anticommutator of any 
two Q charges is some specific polynomial in {J 6 , J 4 , J 2 } whose order grows linearly with g. 
These and the commutators with the even charges are complicated but can be worked out 
in principle. 

6 Outlook 

The model we have considered in this paper is a very special one, and our treatment profited 
from the isometry A 3 ~ D :i . ft will be quite interesting to find out how the structure 
generalizes, firstly, to higher n in A n _i and, secondly, to general Coxeter systems of rank r. 
Clearly, for the angular system three objects are crucial: 

• the ring of Weyl-symmetric polynomials in {£p} and its independent generators 

• maximal abelian subalgebras of the ring of Weyl-symmetric polynomials in { jC ?j } 

• the elementary Weyl-antiinvariant polynomials in {Cij} 

Assuming superintegrability, we conjecture that the number of ring generators equals 2 r—3 
and r —1 for the first and second item, respectively, but one has to go to n =5 (or r= 4) to 
really test the second assertion. Especially useful will be the answer to the third item, as 
different intertwiners can be employed to partially lift the growing degeneracy of the energy 
levels. 

Naturally, an extension to trigonometric, hyperbolic or even elliptic angular Calogero 
systems may be contemplated. Finally, it will be fascinating to introduce a VT deformation 
into the angular system [37, 38], because it will remove the codimension-one singular loci of 
the potential, thereby connect the n\ disjoint particle sectors and also give meaning to g <0 
states, which are needed for the action of the odd Q charges in the spectrum. I 11 fact, it was 
observed that the action of such kind of charges on physical states gets regularized when 
VT deformations are considered, like in the case of the complex trigonometric Poschl-Teller 
potential (which corresponds to the 3-particle case considered here) [39]. 
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A Appendix: 


We collect the more lengthy formulae here. The most complex elementary intertwining 
relation reads 


sy(ff) t(s) 
1V1 6 +6 


_ / 7(9+1) . -p 64 t( 9 + 1 ) 1 -p622 7(9+1) 7(9+1) -p62 7(9+1) , 71,7(9) 

— I ’J f\ "+ 1 66+4 "T" 1 fifi + 2 + 2 "T 1 I 1VJ R 


66 


66^2 


66 


+ ( - 1U4 9+1) 4 9+1) - 364 9+1) 4 9+1) + 168J? +1) 4 g+1) 4 a+1) 

oc 7(9+1) j{g+l) t( 9 + 1 ) 7(9+1) p 36 t( 9 + 1 ) i -p 342 7(9+1) 7(9+1) ~p 3222 7(9+1) 7(9+1) j{g+l) 

OUU 2 +2 J 2 +2 ' 1 66+6 ' 1 66 +4 +2 ' 1 66 +2 +2 +2 


~p 34 7(9+1) 1 ~p 322 7(9+1) 7(9+1) . ^34 j(g+l) , -p 322 7(9+1) 7(9+1) . p 32 r( 9 +l) 1 -p 3 (A 71,7(9) 
1 66+4 "T" 1 fifi +2 J 2 1 ” 1 1 “ 1 fifi T i "T 1 fi« I lv± i 


66 


66^4 


66 


6G°2 


66 


with the abbreviations 


r[* =-(107+ 840), 

Fef = -3(23 + 36#), 

Til = 2(4411 + 13438# + 12432# 2 + 3168c? 3 ), 

= -24(4826 + 33989c? + 71466c? 2 + 64056c? 3 + 25056c/ 4 + 3456c? 5 ), 

= 48(266 + 255# + 36# 2 ), 

T 342 = -48(329 + 306<7 + 60# 2 ), 
r 3222 = 16(271 + 245c? + 48c? 2 ), 

T 34 = 48(4486 + 10765# + 8184c? 2 + 2160c? 3 ), 

T 322 = -16(80 71 + 229 1 9# + 20280# 2 + 5616# 3 ), 

T 32 = 96(23702 + 151661# + 328148# 2 + 317796# 3 + 139680# 4 + 22464# 5 ), 

r 3 ° = -2304(13134+126277#+440735# 2 +759693# 3 +715926# 4 +372744# 5 +99360# 6 +10368# 7 ) . 
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The full expression for the distinguished conserved charge i? 12 = MgM g is 

R\2 = — 12JgJg + 12 (J 6 T 4 + Jg'Ja) ,J‘2 — JQJ2J2J2 + 2J4J4J4 — 1 4 .Z4 J4 J‘2 J2 + 6J4 J 2 J 2 J 2 J 2 
— \ J2J2J2J2J2J2 + + ^ 6 ^ 4 ) + ^it 22 JQJ2J2 + r )^ 2 JaJA J2 + 'Yq 222 JA J2J2J2 


+ 7g 2 " 22 J2J2J2J2J2 + 76^4^2 + 7g 4 </ 4 ^4 + 1^ ZZ 'h-h'^2 + r )^ 22 ^ J2J2J2J2 + 76^6 + 'YfR J a J2 


,,62 


,,44 


422 


,2222 


,,42 


+ 'Yq 2 ^J 2 J 2 J 2 + 7(3 ^4 + Iq 2 J 2 J 2 + 76^2 , 


22 


where the ^-dependent coefficients are 

7e 64 = -4(140 + 39(7 + 36/), 

7 6 22 = Iff - 24(7 + 96/ , 

7 6 442 = 8(202 + 63(7 + 30/), 

7 6 4222 = -§(388 + 65(7 + 48/), 

7 6 22222 = f(31 + 3c7 + 3/), 

7 g 2 = -64(322 + 203 g - 210/ - 216c/ 3 ), 

7 g 4 = 16(553 + 308c/ - 1884/ - 648/), 

7 422 = |(7045 + 4352(7 + 2712/ - 4320/), 

7 2222 = —1(2157 + 944 g + 48c/ 2 - 1296/), 

7 g = -96(3940 + 1469c/ - 8691/ - 4584/ + 5328/ + 3456/), 

7 / = 16(35466 + 18288 g - 97805/ - 54456/ + 56880/ + 44928/), 

7 g2 2 = -f (33074 + 18470/7 - 94827/ - 41064/ + 51984/ + 38016/), 

7 ^ = 192(5864 + 6158<7 - 57481/ - 32799/ + 119796/ + 68112/ - 60480/ - 41472/) 
7 22 = —32(19568+28410c/—165011/—106098/+323232/+146592/—138240/—82944/) 
7 | = -384(448 - 2274 g - 32767/ - 58215/ + 52068/ + 64296/ - 89856/ - 51840/) 
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Finally, the shift, p[ 9 f l) = j^(R[ 9 2 +1) 


— r[ 2 9 ' > ) takes the form 


1 (3+ 1 ) 
4r 12 


= -54 (j; 


(9+1) t(9+1) 

J A 


+ 4 9 + 1 ) 4 9+1) ) - 12 j f 


(9+1) t(s+1) t(s+1) 

o Jo 


6 


+ 294j! 9+1) ji 9+1) Jn S+1) 


_ i 70 7 ( 9 + 1 ) T (g+ 1) r(g+l) t(s+ 1) _ cn t(9+1) t(9+1) 7 ( 9 + 1 ) 7 ( 9 + 1 ) t(9+1) i -.62 t(9+1) t(9+1) 

a 2 J 2 J 2 u aj 2 j 2 j 2 j 2 j 2 I ’ 1 12+6 J 2 


44 t( 9+1 ) t( 9+1) , TV. 

112 +4 +4 ' '/12 +4 


422 7(9+1) j(9+1)j(9+1) ^2222 7(9+1) 7(9+1) 7(9+1) ,(9- 0 ,„(i 7(9+1) 


’2 °2 I Il2 J 2 °2 


Rr’j. 


2 ^2 + '7l2^6 


V7 


1 „42 7(9+1) 7(9+1) ziiz jyy - r^j jyy - r ^) 7191-17 , 4 7191-1; , 7191-1; 7191-1; , 2 7V9tij , u 

~+ 712+4 +2 ~+ 7l2 +2 +2 + 2 ~+ c/12+4 ~+ // 12+2 +2 ~+ //12+2 ~+ c/12 , 


,.222 7 ( 9 + 1 ) 7 ( 9 + 1 ) 7 ( 9 + 1 ) 


(9+1) 


„22 7 ( 9 + 1 ) 7 ( 9 + 1 ) 


7(9+1) 


and the nontrivial coefficients are as follows, 

r / 62 = 36(25 + 304(7 + 144/), 

= -3(1739 + 4008(7 + 1152/), 
rff = 2(737 - 3924 g - 3168c/ 2 ), 
r/ff 2 = 10013 + 23048c/ + 2208/ , 
rf l2 = -24(274 + 8397# + 27252 g 2 + 22032/ + 5184/), 
qf 2 = 12(10096 + 71365 g + 159372/ + 118224/ + 22464/), 
rjff = -4(151692 + 632591# + 717108/ + 116208/ + 19008/), 

rjf 2 = -24(82623 + 628336 g + 1923894/ + 2950776/ + 2270592/ + 767232/ + 124416/), 
qf 2 = 12(1509081 + 9095222# + 18576756/ + 13405872/ + 1368576/ + 442368/ + 82944/), 
qf 2 = -96(2838989+22731865c/+67596615/+87293784/+37690416/-6715008/-1824768/), 
rff 2 = 1728(972216 + 9826815# + 39041714/ + 75013828/ + 65364036/ + 10127664/ 

- 15564096c/ 6 - 4084992/) . 
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i 

(0) (i) 

n n 

(2) 

V i 

(3) 

v\ 

0 

1 c 

c 2 

c 3 

3 

s cs 

c 2 s 

C 3 S 

6 

c 2 -s 2 c(c 2 —3s 2 ) 

c 2 (c 2 —5s 2 ) 

c 3 (c 2 -7s 2 ) 

9 

3c 2 s—s 3 c(c 2 s-s 3 ) 

c 2 (3c 2 s—5s 3 ) 

c 3 (3c 2 s—7s 3 ) 

12 

c 4 —6 c 2 s 2 +s 4 c(c 4 —10c 2 s 2 +5s 4 ) 

c 2 (3c 4 —42c 2 s 2 +35s 4 ) 

c 3 (c 4 —18 c 2 s 2 +21s 4 ) 

15 

5c 4 s—10c 2 s 3 +s 5 c(3c 4 s—10c 2 s 3 +3s 5 ) 

c 2 (3 c 4 s—14c 2 s 3 +7s 5 ) 

c 3 (5c 4 s—30c 2 s 3 +21s 5 ) 

18 

c 6 —15 c 4 s 2 +15c 2 s 4 —s 6 c(c 6 —21c 4 s 2 +35c 2 s 4 —7s 6 ) 

c 2 (c 6 —27 c 4 s 2 +63c 2 s 4 —21s 6 ) 

c 3 (5c 6 —165c 4 s 2 +495c 2 s 4 —231s 6 ) 

21 

7c 6 s—35c 4 s 3 +21c 2 s 5 —s 7 c(c 6 s—7c 4 s 3 +7c 2 s 5 —s 7 ) 

c 2 (5c 6 s—45c 4 s 3 +63c 2 s 5 —15s 7 ) 

c 3 (5c 6 s—55c 4 s 3 +99c 2 s 5 —33s 7 ) 


Tabic 1: Low-lying wave functions vf' > of the Poschl-Teller model. Notation: s = sin(30) and c = cos(30) 





£ 

4 

4 

h (0) 

0 

0 

0 

{000} 

3 

1 

0 

{111} 

4 

0 

1 

{400} - 3{220} 

6 

2 

0 

{600} - 15{420} + 30{222} 

7 

1 

1 

3{511} - 5{331} 

8 

0 

2 

{800} - 28{620} + 35{440} 

9 

3 

0 

9{711} - 63{531} + 70{333} 

10 

2 

1 

{1000} - 45{820} + 42{640} + 504{622} - 630{442} 

11 

1 

2 

5{911} - 60{731} +63{551} 

12 

4 

0 

36{1200} - 2376{1020} + 2445{840} + 46125{822} + 4893{660} - 215250{642} + 179375{444} 

12 

0 

3 

101{1200} - 6666{1020} + 47100{840} + 8685{822} - 42609{660} - 40530{642} + 33775{444} 


£ 

4 

4 

h {1) 

0 

0 

0 

{000} 

3 

1 

0 

{111} 

4 

0 

1 

3{400} - 11{220} 

6 

2 

0 

3{600} - 39{420} + 196{222} 

7 

1 

1 

5{511} - 13{331} 

8 

0 

2 

{800} - 20{620} + 23{440} + 12{422} 

9 

3 

0 

3{711} - 27{531} + 56{333} 

10 

2 

1 

15{1000} - 425{820} + 576{640} + 7568{622} - 14454{442} 

11 

1 

2 

35{911} - 476{731} + 477{551} + 204{533} 

12 

4 

0 

12{1200} - 456{1020} + 657{840} + 13581{822} + 1137{660} - 88842{642} + 114007{444} 

12 

0 

3 

813{1200} - 30894{1020} + 165652{840} + 72131{822} - 147943{660} - 169702{642} + 57527{444} 


£ 

4 

4 

h {2) 

r %A 

0 

0 

0 

{000} 

3 

1 

0 

{111} 

4 

0 

1 

5{400} - 19{220} 

6 

2 

0 

5{600} - 63{420} + 506{222} 

7 

1 

1 

{511} -3{331} 

8 

0 

2 

35{800} - 644{620} + 729{440} + 552{422} 

9 

3 

0 

21{711} - 207{531} + 634{333} 

10 

2 

1 

{1000} - 25{820} + 38{640} + 592{622} - 1374{442} 

11 

1 

2 

21{911} - 300{731} + 303{551} + 200{533} 

12 

4 

0 

980{1200} - 31752{1020} + 51045{840} + 1248957{822} + 84757{660} - 9568482{642} + 16884639{444} 

12 

0 

3 

411005{1200} - 13316562{1020} + 65753340{840} + 36003897{822} - 58071653{660} - 80991402{642} + 7379439{444} 


Table 2: Some wave functions v[ 9) = r e 69 A 9 h < f' > 


of the terahexahedric model, {rst} 
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